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In the framework of various Glauber-like models we compute several 
correlation observables in nuclear collisions at the SPS, RHIC, and LHC 
energies. We analyze fluctuations of the eccentricity of the fireball created 
in the collision, in particular the variable-axes harmonic moments e* , as well 
as the fluctuations of multiplicity of charged particles. We find moderate 
model dependence of the scaled standard deviation a(e*)/e* on the choice 
of the particular Glauber model. For all considered models the values of 
er(e*)/e* range from ~ 0.5 for central collisions to ~ 0.3-0.4 for peripheral 
collisions. The results are confronted to the recent measurement of the 
elliptic- flow fluctuations at RHIC. We also find that the dependence of 
multiplicity fluctuations on the centrality of the collision is too weak to 
explain the measurements at the SPS energies. The magnitude of the 
Glauber multiplicity fluctuations increases by about 20% from the RHIC 
to LHC energies. 

PACS numbers: 25.75.-q 

1. Introduction 

Studies of correlations and fluctutations have become a major part of 
the heavy-ion-collisons physics program, as these observables may carry 
valuable information on the dynamics of the system. While the dynamical 
nature of correlations is of great theoretical interest, a part of the observed 
effect originates from purely statistical phenomena, such as the fluctuation 
of the number of participants in a given centrality class, fluctuations of the 
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shape and orientation of the fireball, etc. It is thus important to understand 
in detail the "background" of these non-dynamical fluctuations. 

In this talk we present predictions of a variety of Glauber-like models of 
the initial stage of a heavy-ion reaction for several correlation observables 
studied at SPS and RHIC, which can also be measured in the forthcom- 
ing LHC experiments. We analyze fluctuations of the eccentricity in the 
transverse plane of the initial fireball, including the variable-axes harmonic 
moments, e*. The predictions are compared to the recent measurements of 
the fluctuations of the elliptic-flow coefficient at RHIC. 

New results are shown for the fluctuations of multiplicity of charged 
particles, and compared to the NA49 results. It is found that the fluctua- 
tion induced by the Glauber-model are not sufficient to explain these data, 
leaving room for dynamical effects. 

The formalism used in this talk is described in detail in Refs. [IJ [2]. In 
particular, all details concerning the statistical methods and the variants of 
the Glauber models may be found there. 

2. Simulations in Glauber-like models 

With the help of a computer program GLISSANDO [2] we have analyzed 
several variants of Glauber-like Monte-Carlo models: 

• The standard wounded nucleon model [3] . The wounding cross section 
a w is equal 32 mb, 42 mb, and 65 mb @] for SPS, RHIC, and LHC 
energies, respectively. 

• The mixed model, amending wounded nucleons with some admixture 
a of binary collisions [3 16] . The successful fits to particle multiplicities 
(see Ref. [BJ) give a = 0.145 at ^/sTn = 200 GeV and a = 0.12 at 
\j $nn = 17.3 GeV. For LHC energy y / s/v r /v = 5.5 TeV we made an 
educated guess for the mixing parameter, a = 0.2. 

• We also analyze a model with hot spots (see Ref. [7]), assuming that 
the cross section for a semi-hard binary collisions producing a hot-spot 
is tiny, <7hot-spot = 0.5 mb for all energies, however when such a rare 
collision occurs it produces on the average a very large amount of the 
transverse energy equal to acr w /<7hot-spot- 

• Each source from the previously described models may deposit the 
transverse energy with a certain probability distribution. Thus, we su- 
perimpose the r distribution, g(w, k), over the distribution of sources 
g(w,K) = w k ~ 1 k k exp(— kw)/T(k). Here we do this superposition on 
the hot-spot model, labeled hot-spot+T . We set k = 0.5, which gives 
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var(u>) = 2 for the analysis of elliptic flow fluctuations. In order to 
reproduce multiplicity distribution in p+p collisions at SPS energies, 
we set k = 1, which gives var(io) = 1 for analysis of multiplicity 
fluctuations. 

The four considered models (wounded-nucleon, mixed, hot-spot, and hot- 
spot+r) differ substantially in the number of sources and the amount of the 
built-in fluctuations. 



3. Event-by event fluctuations of the elliptic flow 

The proper description of the mechanism of the fluctuations of the el- 
liptic flow may be very helpful in understanding the dynamics of heavy-ion 
collision, especially in its early stage [8]. Moreover, for the first time the 
elliptic-flow fluctuations have recently been measured at RHIC [9| flOl \TT\ . 
The interpretation of these results connects the fluctuations of the eccentric- 
ity coefficient, e*, with the fluctuations of the variable-axes flow coefficient, 
denoted in this talk as For sufficiently small elliptic asymmetry, which 
is an experimental fact, one expects on purely hydrodynamic grounds the 
relation 



a(v*) _ a(e*) 
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Fig. 1. (Color online) Fluctuations of e* in two variants of the Glauber models, 
compared to the data from Refs. [9], |T0l [TTJ . See the text for details. 
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Comparison of the data to our Glauber calculations is made in Fig. [TJ For 
central collisions we expect the result following from the central limit theo- 
rem p] for uncor related sources, 

a{vl)/vl ~ a(e*)/e* = ^/A/tt - 1 ~ 0.52 (6 = 0), (2) 

which is compatible with the data, although the error bars are large. We 
note that, with the identification ([T]), the Glauber models reproduce prop- 
erly the data for not-too-peripheral collisions, where the approach is cred- 
ible. Other variants fall between the wounded- nucleon model, which has 
the lowest amount of fluctuation, and the hot-spot+r model, which has 
the strongest fluctuations. More accurate measurements are needed to dis- 
criminate the model predictions. We note, however, that the statistical 
fluctuations built in the Glauber model are sufficient to explain the data. 

In Fig. [2] we compare the results of simulations of the elliptic- flow fluc- 
tuations at SPS and RHIC with the predictions for the LHC energies. The 
calculation is done for two extreme variants of Glauber-like models, namely 
the wounded nucleon model, where the fluctuations of initial shape are 
smallest, and the model with hot-spots+r, where the initial shape fluctu- 
ates are most prominent. For the most central and most peripheral collisions 
the results are very similar for both models at all energies, however at inter- 
mediate centralities there are much higher fluctuations for the LHC energies 
with the hot-spots+r model than in other expected. 




Fig. 2. (Color online) Fluctuations of u| at SPS, RHIC, and LHC energies from the 
Glauber model with wounded nucleons only (full lines) and the hot-spot+L model 
(dotted lines) . 



epiphany 'rybczynski printed on March 29, 2008 



5 



4. Multiplicity fluctuations 

Recently the NA49 Collaboration published results [12] on centrality 
and system-size dependence of multiplicity fluctuations observed in Pb+Pb 
minimum bias and in p+p collisions. Unexpectedly, the scaled variance 
Var(A r )/(A r ), where Var(iV) is the variance and (N) the average multiplic- 
ity of the observed charged particles, changes non-monotonically when the 
number of wounded nucleons grows. The scaled variance is close to unity in 
peripheral and central collisions, however it shows a very prominent peak 
at N w ~ 70. The measurement has been performed at the top SPS collision 
energy ^/snn = 17.3 GeV in the transverse momentum and pion rapidity 
intervals (0.005, 1.5) GeV and (1.1,2.6), respectively. The azimuthal accep- 
tance has been also limited, and only about 20% of all produced negative 
particles have been used in the analysis. In Fig. Owe show NA49 results for 
the negatively-charged particles compared to the Glauber-like simulation 
of the model hot-spots+r in the NA49 azimuthal acceptance. The experi- 
mental acceptance of 20% has been implemented in the simulation. Despite 
the fact that the fluctuations in the hot-spot+r model are the highest of 
all considered Glauber models, it is evident from the plot that the Glauber 
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Fig. 3. (Color online) Scaled variance of the multiplicity distribution of negatively 
charged hadrons produced in p+p and Pb+Pb collisions at the top SPS energy, 
plotted as a function of centrality determined by the number of projectile partic- 
ipants measured by NA49 experiment. The line is obtained with the hot-spot+r 
model with proper implementation of the experimental acceptance. The data are 
from Ref. [T§]. 
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calculation alone cannot describe the non-monotonic dependence of scaled 
variance on centrality. Near the peak the model falls below the data, while 
at larger N w it goes above the data. For other models the obtained curve 
is even lower. In our view dynamical effects must be incorporated in order 
to understand the phenomenon |X3|, HH] . 

In Fig. [J] we compare the centrality dependence of the model multiplicity 
fluctuations for the SPS, RHIC, and LHC energies in the full acceptance. 
Here we use scaled standard deviation of multiplicity distribution as a mea- 
sure of multiplicity fluctuations. Multiplicity fluctuations at SPS and RHIC 
energies are very close to each other, however we notice an increase of for the 
LHC by about 20%. This effect is due to the larger value of the wounding 
cross section a w . The conclusion is qualitative, since as mentioned above, 
the Glauber models do not explain the details of the SPS measurement of 
the multiplicity fluctuations. 

5. Conclusions 

Our main results are as follows: 

• The fluctuations of the eccentricity of the fireball are related to the 
fluctuations of the elliptic flow coefficient. They depend rather moder- 
ately on the chosen Glauber model. Our results agree with the recent 




Fig. 4. (Color online) Scaled standard deviation of the multiplicity distribution 
of negatively charged hadrons plotted as a function of the number of projectile 
participants for SPS, RHIC and LHC energies, all results for the Glauber Monte- 
Carlo hot-spot+r model. 
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PHOBOS and STAR measurements of the V2 fluctuations. The agree- 
ment indicates that this quantity is dominated by the Glauber-model 
statistics. 

• The Glauber-model predictions for the multiplicity fluctuations fail 
short of the experimental results of Ref. |12j . leaving room for strictly 
dynamical effects [HI [14] . Such effects should be introduced in order 
to understand the non-monotonic dependence of the scaled variance 
of multiplicity on the number of produced particles. 
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